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Abstract. With the exception of confinement the three-dimensional Nambu— Jona-Lasinio (NJL3) model 
incorporates many of the essential properties of QCD. We discuss the critical properties of the model at 
nonzero temperature T and/or nonzero chemical potential jj,. We show that the universality class of the 
thermal transition is that of the d — 2 classical spin model with the same symmetry. We provide evidence 
for the existence of a tricritical point in the (/i, T) plane. We also discuss numerical results by Hands et 
al. which showed that the system is critical for /i > /ic and the diquark condensate is zero. 

PACS. 71.10.F Lattice fermion models - 05.70.F Phase transitions 



1 Introduction 

Phase transitions in QCD at nonzero temperature and/or 
nonzero baryon density have been studied intensively over 
the last decade both analytically and numerically. How- 
ever, since the problem of chiral symmetry breaking and 
its restoration is intrinsically non-perturbative, the num- 
ber of available techniques is limited and most of our 
knowledge about the phenomenon comes from lattice sim- 
ulations. Because of the complexity of QCD with dynami- 
cal fermions, studies so far have been done on lattices with 
modest size and in various cases the results are distorted 
by finite size and discretization effects. 

The NJL model has been proved to be an interesting 
and tractable laboratory to study chiral phase transitions 
both numerically by means of lattice simulations and an- 
alytically in the form of large- TV/ expansions [1-9]. The 
lagrangian density of the J7(l)-symmetric model is 



C = tpi{(^ + m + a + ij5n)ipi + 



(1) 



where the index i runs over Nj fermion flavors. There are 
several features which make this model interesting for the 
modelling of strong interactions: (i) The spectrum of ex- 
citations contains both "baryons" and "mesons" , namely 
the elementary fermions / and the composite // states 
|ll0| ; (ii) For sufficiently strong coupling 5^ > g"^ it exhibits 
spontaneous chiral symmetry breaking implying dynam- 
ical generation of a fermion mass m/, the pion field tt 
being the associated Goldstone boson; (iii) For 2 < d < 4 
there is an interacting continuum limit at a critical value 
of the coupling, which for d — 3 has a numerical value 
Ri 1.0 in the large- A^/ limit if a lattice regularisa- 
tion is employed |3] . There is a renormalisation group UV 
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fixed point at — gl., signalled by the renormalisability of 
the 1/Nf expansion entirely analogous to the Wilson- 
Fisher fixed point in scalar field theory; (iv) Numerical 
simulations with baryon chemical potential /i 7^ show 
qualitatively correct behavior, in that the onset of matter 
occurs for u of the same order as the constituent quark 
scale TO/ 1^, rather than for ^ k, to^/2, which happens in 
gauge theory simulations with a real measure det(M'l'M) 
because of the presence of a baryonic pion in the spec- 
trum. This makes NJL3 an ideal arena in which to test 
strongly interacting thermodynamics. 

In section we discuss the universality of the T 7^ 
transition [^||]. In section ^ we present results from a 
study of the phase diagram in the (/i, T) plane [§ which 
support the existence of a tricritical point on the line that 
separates the chirally broken from the chirally symmetric 
phases. We also discuss numerical results which support 
the non-existence of a diquark condensate for /i > /ic @ . 



2 Universality at nonzero temperature 

Although there is little disagreement that the chiral phase 
transition in QCD with two massless fermions is second 
order no quantitative work or simulations have been done 
that decisively determine its universality class. Universal- 
ity arguments are very appealing due to their beauty and 
simplicity. In essence they can be phrased as follows: At 
finite T phase transitions, the correlation length diverges 
in the transition region and the long range behaviour is 
that of the {d—1) classical spin model with the same sym- 
metry, because the IR region of the system is dominated 
by the zero mode of the bosonic field and the contribution 
of non-zero modes does not affect the critical singulari- 
ties but can be absorbed into non-critical, non-universal 
aspects of the transition. Consequently, fermions, which 
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Table 1. Summary of FSS results and comparison with 2d 
Ising and mean field scaling behaviour. 



Exponents 


FSS 


Z(2) 


MP 


V 


1.00(3) 


1 


0.5 


13m./ y 


0.12(6) 


0.125 


1 


-y/v 


1.66(9) 


1.75 


2 



satisfy anti-periodic boundary conditions, and do not have 
zero modes are expected to decouple from the scalar sec- 
tor. In another language, the classical thermal fluctuations 
whose energy is 0{KbT) dominate over quantum fluctua- 
tions with energy 0{Tujj) for soft modes of the field and the 
effective {d — l)-dimensional theory for the bosonic flelds 
near Tc is a classical statistical theory. A possible loophole 
to this standard scenario is that the mesons are composite 
// states and their size and density increase as T ^ Tc. 
Therefore, if the transition region can be described as a 
system of highly overlapping composites the violation of 
the bosonic character of the mesons may be maximal and 
the fermions become essential degrees of freedom irrespec- 
tive of how heavy they are. At leading order in 1/Nf the 
model has a second order phase transition at Tc — -^^^2 
1^ with Landau- Ginzburg mean field scaling. We studied 
the critical behavior of the Z2-symmetric model with fi- 
nite Nf by performing lattice simulations in the vicinity 
of the critical point. The temporal lattice size was Lt = 6 
and the spatial size varied from Ls = 18 to 50. The expec- 
tation value of the auxiliary signia field {S = {a)) serves 
as a convenient order parameter for the theory's critical 
point. We simulated the model exactly for Nf — 12 with 
the Hybrid Monte Carlo method. The staggered fermion 
lattice action and further details concerning the algorithm 
can be found in By using the finite size scahng (FSS) 
method we extracted the critical exponents. The results 
which are summarized in Table 1 support the dimensional 
reduction scenario, because the values of the exponents 
are in good agreement with those of the 2d Ising model 
rather than the mean field theory ones. 

Next we tried to understand how the large- A^/ mean 
field theory prediction reconciles with the dimensional re- 
duction and universality arguments. The answer is that 
the large-A^y description has its applicability region. As 
we discuss in detail for a Yukawa theory in |^ the phe- 
nomenon which leads to an apparent contradiction is the 
suppression of the width of the non-mean field critical re- 
gion by a power of 1/Nf. The 2d Ising critical behavior 
sets in when T ^ ma{T) {ma{T) is the thermal mass of 
the (7 meson). If the cutoff yl ^ T the renormalized self- 
interaction coupling A(T) in the large- A^/ limit is close to 
the IR fixed point of the Yukawa theory and is given by 
A(T) ~ T^^'^/Nf for 2<d<4. The mean field approxima- 
tion breaks down because of self-inconsistency when the 
value of the coupling of the {d—1) dimensional scalar the- 
ory \d-i = TX{T) on the scale m^~''(T) (the power d - 5 
comes from dimensional analysis) is not small anymore. 
Therefore, for d — 3 the Ginzburg criterion for the ap- 
plicability of the mean field scaling is given by m^lT) ^ 
T / y/Nf. This scenario was verified numerically in 
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Fig. 1. E{T,ii)/E{T,0) vs. fi/nif at different values of T. 
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Fig. 2. ri{fi)/mf vs. i-i/mf at different values of T. 

Additional evidence in favor of the dimensional re- 
duction scenario was produced in studies of the U{1)- 
symmetric NJL3 model Both analytical and numer- 
ical results showed that its phase structure at T ^ is 
the same as the 2c? XY model. It has two different chi- 
rally symmetric phases, one critical and one with finite 
correlation length, separated by a Berezinskii-Kosterlitz- 
Thoulcss transition. 



3 Results at non-zero chemical potential 

The action of the NJL model remains real even after the 
introduction of nonzero chemical potential /i, which means 
we can study the physics of the high density regime us- 
ing standard Monte Carlo techniques. In the presence of 
a Fermi surface with Fermi momentum pp the creation 
of // pairs with zero net momentum is suppressed, be- 
cause one can only create particles with p > pp. So as the 
fermion number density 7y(/i) grows the chiral symmetry 
breaking is suppressed. The large- A^/ description of the 
fj. ^ chiral phase transition predicts a first order transi- 
tion for T = and a continuous transition for T > [Q. 
Furthermore, the critical value of the chemical potential 
fj,c is equal to the value of the fermion mass at ^ = 0, 
which indicates that materialization of the fermion itself 
drives the symmetry restoration transition. Interactions as 
expected decrease fic below the mean field result [|. Work 
by Stephanov suggests that any non-zero density sim- 
ulation which incorporates a real path integral measure 
proportional to det (MM^) is doomed to failure due to the 
formation of a light baryonic pion from a quark q and a 
conjugate quark q'^. The NJL model however, does not ex- 
hibit such a pathology, because the realization of the Gold- 
stone mechanism in this model is fundamentally different 
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Fig. 3. Fermion dispersion relation at /i = 0.8. 



from that in QCD. In the NJL model the Goldstone mech- 
anism is realized by a pseudoscalar channel pole formed 
from disconnected diagrams and the connected diagrams 
yield a bound state of mass « 2m/. This implies the ab- 
sence of a light qq'^ state. As expected our simulations 
of the Z2 symmetric NJL3 with Nf = 4 |^ did not pro- 
vide any evidence for the existence of a nuclear liquid-gas 
transition at /i < /ic- It was shown in | [Tl[ | that in the 
NJL model there is no saturation density for stable mat- 
ter. In order to get the saturation features the authors of 
lini introduced a chirally invariant scalar-vector interac- 
tion term which cures the binding problem. Furthermore, 
our results showed that the second order nature of the 
r ^ 0, /.t = transition remains stable down to low T and 
large /i. In Fig.Q we plot the normalized order parameter 
S{T, fj)/ S{T,Q) as a function of n/ruf at different val- 
ues of T. It is clear from the shapes of these curves that 
the transition becomes sharper as we decrease the tem- 
perature. In Fig.^jwe plot the normalized fermion number 
density as a function of the chemical potential at different 
values of T. In the limit T ^ we see that the fermion 
density is strongly suppressed before the transition and 
then jumps discontinuously. By performing detailed finite 
size scaling analysis which allowed us to distinguish be- 
tween second order and weak first order transitions we 
showed that the tricritical point lies on the section of the 
phase boundary defined by T/T^ < 0.23, ^l/^lc > 0.97 g. 
This result shows that higher order corrections in the 1 /Ay- 
expansion for the nature and the location of the transition 
points in the phase diagram are small in this model. 

It is well-known that at high density the diquark con- 
densate is nonzero in models of strongly interacting matter 
which either assume that the interaction between quarks 
is due to one-gluon exchange Q, or by using effective 
four-fermion vertices resulting from the presence of instan- 
tons in the QCD vacuum [ |l3[ . Unfortunately, theoretical 
studies of color superconductivity are limited to perturba- 
tive and self-consistent methods, because of the notorious 
"sign problem" in QCD. Hands et al. studied numerically 
the SU{2) (g) S't/(2)-symmetric NJL3 and found no evi- 
dence for a condensate (qq) 7^ in the model's high den- 
sity phase. Their results with a nonzero diquark source are 
consistent with a critical behavior {qq{j)) oc through- 
out the dense phase with S falling in the range 3 — 5 for the 
/z values studied. This suggests that the model is a two- 
dimensional superfluid as first described by Kosterlitz and 



Thouless for thin films of "^He but with the universality 
class determined by the presence of massless relativistic 
fermions. Results for the dispersion relation E{k) in the 
spin-^ sector are shown in Fig.^. For k < kp the lowest 
excitations vacate states in the Fermi sea, and hence are 
'hole-like', whereas for k > kp, excitations add quarks to 
the system and are 'particle-like'. There is no sign for any 
discontinuity on the Fermi surface characteristic of a BCS 
gap Z\ ^ in NJL3. However, recent numerical results 
provided evidence that {qq} 7^ in NJL4 |14]. 



4 Summary 

We discussed the basic features of the phase diagram of 
NJL3. The universality class of the T 7^ transition is that 
of the d = 2 classical spin model with the same symmetry. 
The non-trivial critical region of the Z2-symmetric model 
is suppressed by a factor 1 / y^Nf. The non-zero density 
phase transition is strongly first order and the lattice sim- 
ulations provided evidence for the existence of a tricritical 
point on the critical line at small T and large fj,. The sim- 
ulations showed no evidence for the existence of non-zero 
diquark condensate at T = and /j, > fj,c. The results are 
consistent with a critical behaviour throughout the dense 
phase, suggesting that the model is a two-dimensional rel- 
ativistic superfluid. 
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